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Abstract 

Atomic Compton profiles (CPs) are a very important property which provide us information 
about the momentum distribution of atomic electrons. Therefore, for CPs of heavy atoms, rela- 
tivistic effects are expected to be important, warranting a relativistic treatment of the problem. In 
this paper, we present an efficient approach aimed at ab initio calculations of atomic CPs within 
a Dirac-Hartree-Fock (DHF) formalism, employing kinetically-balanced Gaussian basis functions. 
The approach is used to compute the CPs of noble gases ranging from He to Rn, and the results have 
been compared to the experimental and other theoretical data, wherever possible. The influence of 
the quality of the basis set on the calculated CPs has also been systematically investigated. 
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I. INTRODUCTION 



Recent years have seen tremendous amount of progress in the field of relativistic electronic 
structure calculations of atoms and molecules using Dirac-equation-based approaches^. 
Particularly noteworthy are the advances made in the field of basis-set-based relativistic 
electronic structure theory pioneered by Kim^, and Kagawas}. Although, initially, the 



basis-sets employed in the calculations were of the ordinary Slater- type 2, 3], however, now- 
a-days, the preferred basis functions are those which incorporate the so-called kinetic-balance 



aa. 



condition between the large and the small component basis functions [4|, 15|, 16||. The most com- 
monly used variety of such functions in relativistic electronic-structure calculations are the 
kinetically-balanced Gaussian functions (KBGFs) which have not only been instrumental in 
avoiding the problem of 'variational collapse', but have also allowed the import of efficient 
algorithms developed in basis-set-based nonrelativistic quantum chemistry. Using such ba- 
sis functions, calculations are now routinely performed both at the mean-field Hartree-Fock 



a a, 



as well as at the correlated leve 



em- 



(henceforth Dirac-Hartree-Fock (DHF)) level 
ploying methods such as the configuration-interaction (CI) approach, both for atoms[9|, and 
molecules flo| . 

However, the progress in calculating wave functions and atomic energies using KGBFs 
has not been matched by the progress in computing expectation values corresponding to 
various physical quantities. For example, atomic Compton profiles (CPs) are a very im- 
portant property which provide us information about the momentum distribution of atomic 
electrons, and help us in interpre ting the x-ray Compton scattering data from atoms in the 



large momentum-transfer regime 



Compton profiles are also very useful in understand- 



ing the bonding properties, as one makes a transition from the atomic scale to the scale of 



condensed matter 



Ufl. Indeed, the nonrelativistic Schrodinger equation based calculations 



of CPs of atomic and molecular systems both within an ab initio, as well as model-potential 



based, formalisms are quite well developed 



As recently demonstrated by us, and several 



other authors earlier on, that such nonrelativistic ab initio calculations of CPs can also be 



performed on crystalline systems 



121 ] . However, for systems involving heavy atoms, on intu- 



itive grounds one expects that the relativistic effects will become quite important, the reby 



requiring a relativistic treatment of the problem Long time back Mendelsohn rf at. 



and Bigss et a/.[15y presented the first fully-relativistic calculations of atomic CPs which 
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were performed at the DHF level, employing a finite-difference based numerical approach. 
Yet, since that time, there has been hardly any activity in the field, which is surprising 
given the fact that now relativistic electronic structure calculations are routinely performed 
employing KBGF basis functions. Therefore, in this work, our aim is to report the first 
calculations of atomic CPs at the DHF level, employing a basis set composed of KBGFs. 
Our approach is based upon analytic formulas for the CP matrix elements with respect to 
a KBGF basis set, whose derivation is presented in the Appendix. The DHF calculations 
of atomic CPs are presented for the entire rare gas series (He to Rn), and our results are 
compared to experimental data, wherever available. Additionally, our results for Ar, Kr, Xe, 
and Rn are also compared to the DHF results of Mendelsohn et a/. jl^J], and Bigss et al. 15], 
and excellent agreement is obtained between the two sets of calculations. 

At this point we would like to clarify one important aspect related to the relativistic 
effects which our calculations are computing, in light of the fact that there have been several 
papers in the literature dealing with a relativistic treatment of Compton scattering of bound 

nnn 

electrons [16|, [T7J, H8fl. Several authors have pointed out that for very large photon energies, 



a fully relativistic treatment, within the framework of quantum-electrodynamics, of the 



Compton scattering from bound electrons is essential 



13j. When such a treatment of the 



problem is performed, it is not clear whether the Compton scattering cross-sections can at all 



be written in terms of Compton profiles 



17 



18fl. Our work presented here, however, 



does not correspond to that regime of photon energies. What we mean by the relativistic 
effects here are the changes in the computed CPs because of a relativistic treatment of 
the bound electrons within a Dirac Hamiltonian based formalism. Thus, our calculations 
assume that the Compton scattering from atomic electrons can be described in terms of 



the CPs under the impulse approximation 



19(| . The electron momentum densities needed to 



calculate the CPs, however, are computed from the Dirac orbitals of the atomic electrons. 
This approach is identical to the one adopted in the earlier DHF calculations \\\ 

The remainder of this paper is organized as follows. In section [Til we present the basic 
theoretical formalism behind the present set of calculations. Next in section [TTT] we present 
and discuss the results of our calculations. Finally, in section [TV] our conclusions, as well 
as possible future directions for further work are discussed. Additionally, in the Appendix 
we present the derivation of the closed-form formulas for CPs over KBGFs, used in our 
calculations. 
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II. THEORY 



Our theory is based upon the Dirac-Coulomb Hamiltonian 

H = ]T(ca 4 • p, + c 2 (A - 1) + Vnudn)) + £ — , (1) 

i i<j l l 

where c is the speed of light, p is the momentum operator, Vnuc( r ) is the electron-nucleus 
interaction potential, indices i and j label the electrons of the atom, and is the distance 
between the ith and jth electrons. For Vnuc(r) a spherical finite-nucleus approximation 



is employed, with the radius estimated as 2.2677 x 10 5 A 1 / 3 , where A is the atomic mass 



rfl- 



/ 



number |7j. The Dirac matrices are chosen to be a = and (3 




( 



where 0, 



V 



i o 

o -I 

I, and a, represent the 2x2 null, identity and Pauli matrices, respectively. Eq. (DD) is solved 
under the DHF approximation utilizing spherical symmetry with the orbitals of the form 

/ 



-1 



(2) 



P nK (r)XKm(0,(j)) 
\iQn K {r)x-Km(Q,<f>) 

where P nK {r) and Q nK (r) are the radial large and small components, and Xnm(9, (ft) is the two- 
component angular part composed of Clebsch-Gordon coefficients and spherical harmonics. 
In the basis-set approach adopted here, the radial parts of the wave function are expressed 
as linear combination of radial Gaussian type of functions 

^«(r) = E c £aS(0, 

i 

and 

Q nK (r)=j:C s Ki g s Jr), 

i 

where and C^, are the expansion coefficients of the large and small component basis 
functions, respectively. The large-component basis function is given by 

g L Kl {r) = N L Kl r^e-^\ (3) 

while the small-component basis function is obtained by the kinetic-balancing condition 

9fL = Ni{^ + -)A{r). (4) 
ar r 

Above n K is the principal quantum number associated with a symmetry species {n K = 
1, 2, 2, 3, 3, . . ., for symmetry species s, Pi/%, P3/2, $3/2, ^5/2, - - -), o» is the Gaussian exponent 



of the zth basis function, and iV£, iV£ are the normalization coefficients associated with the 



Kl 1 Kl 



large and the small component basis functions, respectively. 

Under the impulse approximation the differential cross-section of Compton scattering 
of x-rays from many-electron systems is proportional to the Compton profile 

J(<l)=f J dp x dp y p(p), (5) 

where p(p) is the momentum distribution of the electrons before scattering and q is the 
component of the momentum of the electron along the scattering vector, assumed to be 
along the z direction. Under the mean-field DHF approximation, for a closed-shell atom, 
the expression for the CP reduces to 

J(q)=J2(2ji + l)Jn iK M, (6) 

where ji is the total angular momentum of the ith orbital while, J„ iKi is the CP associated 
with it 

J ™b) = 2j q {\ P ^(P)\ 2 + \Qn^(p)\ 2 }pdp, (7) 

where P Hl Ki (p) and Q ni Ki ip) are the Fourier transforms of the radial parts of the large and 
small components, respectively, of the zth occupied orbital (cf. Eq. (j2J)) and are defined as 

PuiKiip) = ( 27r ) 3 /2 J rP " lKi; (r)ji A (pr)dr, (8) 

and 

4tt f°° 

Qu^Xp) = ( 2?r )3/2 J r QniKi(r)ji B (pr)dr, (9) 

where ji A {pr)(ji B (pr)) is the spherical Bessel function corresponding to the orbital angular 
momentum Ia(Ib) of the large (small) component. Therefore, calculation of atomic CPs 
involves computation of two types of integrals: (i) radial Fourier Transforms of Eqs. ([HI) and 
Q, and (ii) momentum integrals of the Fourier transformed orbitals in Eq. (jjj). When one 
solves the DHF equation for atoms using the finite-difference techniques, then, obviously 
the calculation of atomic CPs mandates that both these types integrals be computed by 
means of numerical quadrature. However, for the basis-set-based approach adopted here, 
in order to facilitate rapid computation of atomic CPs, it is desirable to obtain closed-form 
expressions for both types of integrals with respect to the chosen basis functions. Indeed, 
we have managed to derive closed-form expressions for the atomic CPs with respect to the 
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KBGFs, which can be easily computer implemented. It is easy to see that within a KBGF 
based approach, the integral of Eq. ((Jj), can be computed in terms of the following two types 
of integrals 

Jij' K i.(l) = 2j q P9Ki(p)9Kj{p)dp, 

and 

<j 1 f°° 

J ij K ^) = 2 J P9Ki(p)9Kj(p)dp, 

where g^ip) and g^ip) are the radial Fourier Transforms (cf. Eqs. (JSj) and (j2j)) of the large 
and small component basis functions g^(r), and g^(r), respectively. Obtaining closed- form 
expressions for J^' K (q) and J^' K (g) expressions was not an easy task, and those formulas, 
along with their derivation, are presented in the Appendix. Additionally, elsewhere we have 
described a Fortran 90 computer program developed by us, which uses these expressions to 



compute t 
of KBGFs 



re atomic CPs from a set of given Dirac orbitals expressed as a linear combination 



3- 

Here we would like to comment on possible quantitative manifestations of relativistic 
effects in Compton profiles. One obvious way to quantify the relativistic effects on the CPs 
is by comparing the values obtained from the DHF calculations with those obtained from 
nonrelativistic HF calculations. There is another way by which one can judge the influence 
of relativistic effects on Compton profiles, that is by comparing the orbital CPs of different 
fine structure components. For example, in nonrelativistic calculations, np, nd,. . . orbitals 
have only one set of values each for the orbital CPs. However, in relativistic calculations, 
each such orbital splits into two fine-structure components, i.e., npi/ 2 /np 3 /2, ^3/2/^5/2, 
which, if the relativistic effects are strong, can differ from each other in a significant manner. 
Thus, one expects, that under such situations, the orbital profiles of the two fine-structure 
components will also be significantly different. Therefore, we will also examine this "fine- 
structure splitting" of the orbital CPs of various atoms to quantify the relativistic effects. 

III. CALCULATIONS AND RESULTS 

In this section we present our DHF results on the atomic profiles of the rare gases. The 
DHF orbitals of various atoms were computed using the KBGF based REATOM code of 

n 

Mohanty and Clementi[21]. During the DHF calculations the value of the speed of light 
used was c = 137.037 a.u. Additionally, for obtaining the radius of the nucleus for the 



finite- nucleus approximation description of Vnuc( r ) ? values of atomic mass A were taken to 
be 4.026, 20.18, 39.948, 83.80, 131.3, and 222.0 for He, Ne, Ar, Kr, Xe, and Rn, respectively. 
Using the orbitals obtained from the DHF calculations, the atomic CPs were computed using 

n 

our computer program COMPTON[20j|. Next we present our results for the rare gas atoms, 
one-by-one. In order to investigate the basis-set dependence of the CPs, for each atom, two 
types of basis sets were used: (i) a large universal basis set proposed by Malli et al. 22], 
and (ii) a smaller basis set tailor-made for the individual atom. 



He 



For He, DHF calculations were performed with: (i) well-tempered basis set of Matsuoka 



23j employing 12s functions[23j, and (ii) the universal basis set using 22s 



and Huzinaga 

functions [22J. The computed CPs are plotted in Fig. IIII Al as a function of the momentum 
transfer q. The results of our calculations for some selected values of q are presented in 
table [B For the sake of comparison, the same table also contains the nonrelativistic HF 
results of Clementi and Roetti 24], as well as the experimental results of Eisenberger and 



Reed|25l|. Upon inspection of the table, following trends emerge: (i) Our relativistic CPs 
computed with the well-tempered and the universal basis sets are in excellent agreement 
with each other. This implies that the smaller well-tempered basis set is virtually complete, 
as far as the CPs are concerned, (ii) Our DHF CPs are in excellent agreement with the 
nonrelativistic HF CPs of Clementi and Roetti 2^|. This, obviously, is a consequence of the 
fact that the relativistic effects are negligible for a light atom such as He. (iii) Generally, 
the agreement between the theoretical and the experimental CPs is excellent, implying that 
the electron-correlation effects do not make a significant contribution in this case. 



B. Ne 



DHF calculations wereoerfonned for Ne using: (i) (14., 14p) well-tempered basis set „f 
Matsuoka and Huzinaga[23], and the (ii) large (32s, 29p) universal basis set of Malli et al. [22]. 
In order to facilitate direct comparison with the experiments, the valence CPs (excluding the 
contribution from the Is core orbital) obtained from our calculations are presented in tab 



They are also compared to the nonrelativistic HF results of Clementi and Roetti 



24\ 
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Table I: Relativistic (DHF) Compton profi 
compared to the nonrelativistic HF results 



es of He atom computed using various basis functions, 



241 ] . and the experiments 



2511 



g(a.u.) 


J(g)(WT) a 


J(q)(\Jm) b 


J(g)(HF) c 


J(q)(Exp.) d 


0.0 


1.0704 


1.0704 


1.0705 


1.071 ± 1.5% 


0.1 


1.0567 


1.0567 


1.0568 


1.058 


0.2 


1.0171 


1.0171 


1.017 


1.019 


0.3 


0.9557 


0.9557 


0.955 


0.958 


0.4 


0.8782 


0.8782 


0.878 


0.881 


0.5 


0.7910 


0.7910 


0.791 


0.795 


0.6 


0.7003 


0.7004 


0.700 


0.705 


0.7 


0.6111 


0.6112 


0.611 


0.616 


0.8 


0.5270 


0.5270 


0.527 


0.533 ± 2.3% 


0.9 


0.4503 


0.4503 


0.450 


0.456 


1.0 


0.3820 


0.3820 


0.382 


0.388 


1.2 


0.2712 


0.2712 


0.271 


0.274 


1.4 


0.1910 


0.1910 


0.190 


0.188 


1.6 


0.1344 


0.1345 


0.134 


0.129 


1.8 


0.0952 


0.0952 


0.095 


0.092 


2.0 


0.0678 


0.0678 


0.068 


0.069 


2.5 


0.0307 


0.0307 


0.031 


0.030 ± 15% 


3.0 


0.0148 


0.0148 


0.015 


0.013 


5.0 


0.0014 


0.0014 






8.0 


0.0001 


0.0001 






10.0 


0.00003 


0.00003 







a Our DHF results computed using the well-tempered basis set 



''Our DHF results computed using the universal basis set 



c Nonrelativistic HF results from Ref. 
^Experimental results from Ref. 



25] 



2J 
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12 3 4 

q (a.u.) 



Figure 1: DHF Compt 
the universal basis set 



on profiles of He, J(q), computed using the well-tempered basis set[23|, and 



221 ]. as a function of the momentum transfer q. Profiles obtained using the 



two basis functions are virtually indistinguishable. 



classic experiment of Eisenberger 26|], and more recent experiment of Lahmam-Bennani et 



al. 27fl. Additionally, the total Compton profiles of Ne (including the contribution of the Is 
orbital), computed using both the aforesaid basis sets, are plotted in Fig. El 

Upon inspecting table [Til we notice the following trends: (i) profiles computed using two 
different sets are again in very good agreement with each other, implying that both the basis 
sets are essentially complete, fii) our relativistic profiles are in quite good agreement with 



the nonrelativistic HF profiles 24] essentially implying that even in Ne, the relativistic effects 
are quite negligible. As far as comparison with the experiments is concerned, for smaller 
values of q there is slight disagreement with the theory which progressively disappears as 
one approaches the large momentum-transfer regime. This suggests that electron-correlation 
effects possibly play an important role in the small momentum transfer regime. 

Finally we examine the individual orbital CPs of the Ne atom in Fig. [31 The maximum 
contribution to the total CP for small values of momentum transfer comes from the 2s 
orbital, while in the same region, the smallest contribution comes from the Is core orbital. 
The orbital CP of the 2s orbital varies rapidly with respect to q and becomes quite small for 
q > 2 a.u. On the other hand the orbital profile of the Is orbital shows the least dispersion 
with respect to q, and has the largest magnitude in the large q region, as compared to other 

9 



Table II: Relativistic (DHF) valence Compton profiles of Ne atom computed using various basis 
functions, compared to the nonrelativistic HF results j^]. and the experiments 25]. 



g(a.u.) 


J( 9 )(WT) a 


J(q)(\Jm) b 


J{q) (HF) C 


J(q)(E W .) d 


J(g)(Exp.) e 


0.0 


2.5439 


2.5452 


2.548 


2.582 


2.602 


0.1 


2.5363 


2.5375 


2.540 


2.574 


2.593 


0.2 


2.5128 


2.5140 


2.515 


2.558 


2.560 


0.3 


2.4722 


2.4731 


2.475 


2.519 


2.506 


0.4 


2.4129 


2.4133 


2.418 


2.451 


2.435 


0.5 


2.3342 


2.3339 


2.335 


2.359 


2.340 


0.6 


2.2367 


2.2357 


2.236 


2.249 


2.235 


0.7 


2.1224 


2.1210 


2.120 


2.124 


2.099 


0.8 


1.9947 


1.9933 


1.990 


1.986 


1.966 


0.9 


1.8579 


1.8568 


1.855 


1.839 


1.826 


1.0 


1.7166 


1.7159 


1.715 


1.685 


1.690 


1.2 


1.4360 


1.4361 


1.435 


1.394 


1.417 


1.4 


1.1776 


1.1780 


1.171 


1.140 


1.171 


1.6 


0.9533 


0.9537 


0.951 


0.921 


0.975 


1.8 


0.7663 


0.7665 


0.766 


0.749 




2.0 


0.6142 


0.6144 


0.619 


0.608 




2.5 


0.3559 


0.3558 


0.355 


0.355 




3.0 


0.2125 


0.2123 


0.212 


0.225 




3.5 


0.1318 


0.1319 


0.132 


0.156 




4.0 


0.0852 


0.0853 


0.085 


0.102 




5.0 


0.0397 


0.0397 


0.040 


0rG41 





a Our DHF results computed using the well-tempered basis set [2 



''Our DHF results computed using the universal basis set 



c Nonrelativistic HF results from Ref. 
^Experimental results from Ref. 
e Experimental results from Ref. 



26] 
27] 
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Figure 2: DHF Compton profiles of Ne, J(q), computed using the well-tempered basis set 23], and 



the universal basis set 



22J, as a function of the momentum transfer q. Profiles obtained using the 



two basis sets are virtually indistinguishable. All numbers are in atomic units. 

orbital profiles. The behavior of the 2pi/ 2 /2p 3 / 2 orbital profiles is intermediate as compared 
to the two extremes of Is and 2s profiles. These profiles have lesser magnitude compared to 
the 2s profile for q 0, while they vary more rapidly with respect to q, when compared to 
the Is profile. Another pointer to the insignificance of the relativistic effects for Ne is the 
fact that the difference in the values of the 2py 2 and 2p 3 / 2 is quite small for all values of q. 



C. Ar 



Next, we discuss our calculated Compton profiles of Ar. The DHF calculations on Ar 
atom were performed using the follo wing two basis sets: (i) smaller (16s, 16p) well-tempered 
basis set of Matsuoka and Huzinaga 23j, and the (ii) large (32s, 29p) universal basis set of 
Malli et al. 



221] . Calculated total CPs of Ar, for a selected number of q values in the range 
a.u.< q < 15 a.u., are presented in table IHIl The same table also contains the nonrelativistic 
HF results of Clementi and Roetti numerical-orbital-based DHF results of Mendelsohn 
et al. [14], and the experimental results of Eisenberger and Reed 



25] 



Additionally, in Figs. [4] and [5l respectively, we present our total and orbital CPs of 



11 



Table III: Our relativistic (DHF) total Compton profiles of Ar atom computed using various basis 
sets, compared to the relativistic results of other authors Mj], the nonrelativistic HF results 

n 

the experiments[25]. 



2J|, and 



q{a.u.) 


J(q){WT) a 


J(q)(XJm) b 


J(q)(DRF) c 


J(q){RF) d 


J(g)(Exp.) e 


0.0 


5.0471 


5.0543 


5.05 


5.052 


5.058 ± 0.7% 


0.1 


5.0229 


5.0302 


5.03 


5.028 


5.022 


0.2 


4.9473 


4.9539 


4.95 


4.950 


4.917 


0.3 


4.8130 


4.8171 




4.812 


4.749 


0.4 


4.6143 


4.6144 


4.61 


4.608 


4.526 


0.5 


4.3528 


4.3487 




4.369 


4.259 


0.6 


4.0395 


4.0324 


4.03 


4.028 


3.960 


0.7 


3.6928 


3.6854 




3.690 


3.643 


0.8 


3.3343 


3.3288 




3.328 


3.319 


0.9 


2.9842 


2.9814 




2.982 


3.000 


1.0 


2.6576 


2.6573 


2.66 


2.658 


2.697 ± 1% 


1.2 


2.1071 


2.1088 




2.108 


2.164 


1.4 


1.7011 


1.7022 




1.701 


1.753 


1.6 


1.4163 


1.4166 




1.417 


1.461 


1.8 


1.2198 


1.2197 




1.221 


1.264 


2.0 


1.0825 


1.0824 


1.08 


1.084 


1.129 


2.5 


0.8728 


0.8727 




0.873 


0.904 


3.0 


0.7360 


0.7360 




0.736 


0.744 


3.5 


0.6216 


0.6217 




0.621 


0.634 


4.0 


0.5207 


0.5208 


0.521 


0.520 


0.534 ± 2.5% 


7.0 


0.1773 


0.1774 




0.177 


0.181 


8.0 


0.1300 


0.1300 




0.130 


0.137 


9.0 


0.0981 


0.0981 




0.098 


0.104 


10.0 


0.0758 


0.0757 


0.076 


0.075 


0.078 ± 10% 


15.0 


0.0254 


0.0254 




©7625 


0.025 



a our DHF results computed using the well-tempered basis set [2 



''our DHF results computed using the universal basis set 



14 



C DHF results of Mendelsohn et al. 
^Nonrelativistic HF results from Ref. 



220 



based urJdh finite- difference calculations 



241 




4 6 

q (a.u.) 



10 



Figure 3: Orbital Compton profiles of Ne for 2s (solid line), 2p 3 / 2 /^Pi/2 (dashed line), and Is 
(dotted line), plotted with respect to q. Compton profiles of 2p 1 / 2 and 2p 3 / 2 orbitals are virtually 
indistinguishable. These profiles were computed using the universal basis set 



Ar plotted as a function of the momentum transfer q. From Ar onwards, CP results of 
Mendelsohn et al. flj]] exist, which were computed from the DHF orbitals obtained from 
finite-difference-based calculations. If our calculated CPs are correct, they should be in good 

n 

agreements with those of Mendelsohn et a/.[lj]. Therefore, it is indeed heartening for us to 



note that our CP results computed with the universal basis set 22m are in perfect agreement 

n 

with those of Mendelsohn et a/.[14j] to the decimal places, and for the q points, reported 
by them. As a matter of fact even our CPs obtained using the smaller well-tempered basis 



set 



23| . disagree with those of Mendelsohn et a/.[l4j by very small amounts. Thus, this gives 



us confidence about the essential correctness of our approach. 

When compared to the experiments, for q = 0, our value of CP of 5.054 computed with 



25fl. For 



universal basis set, is in excellent agreement with the experimental value of 5.058 
0.1a.u.< q <0.8a.u. our results begin to overestimate the experimental ones slightly. For 
q >0.9a.u., however, our theoretical results underestimate the experimental results by small 



amounts. The nonrelativistic HF resultsl2j| also exhibit the same pattern with respect to the 

n 

experimental results. Upon comparing our CPs to the nonrelativistic HF CPs[2j|, we notice 
that the two sets of values differ slightly for smaller values of q. However, the difference 
between the two begins to become insignificant as we approach larger values of q, suggesting 
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5 10 

q (a.u.) 



15 



Figure 4: DHF Compton profiles of Ar, J(q), computed using the well-tempered basis set 23], and 
the universal basis set {22], as a function of the momentum transfer q. Profiles obtained using the 
two basis sets can be seen to differ slightly for q m 0. 

that the relativistic effects will be most prominent for q ~ 0. 

Finally we examine the contributions of the individual orbitals to the atomic CP in Fig. 
El which presents the orbital profiles of all the orbitals of Ar. We observe the following 
trends: (i) 3s profile has the maximum value at q — 0, followed by 3^3/2/3^1/2 profiles. 
The minimum value at q = corresponds to the Is profile, (ii) Profiles of outer orbitals 
vary more rapidly with q, as compared to the inner ones. In other words, profile flattening 
occurs as one moves inwards from the valence to the core orbitals. (iii) Again no significant 
fine-structure splitting is observed, in that the profiles of np 3 / 2 and npi/2 orbitals differed 
from each other by small amounts, pointing to the smallness of relativistic effects. 

D. Kr 

Now, we discuss our DHF results of Compton profile of Kr. The DHF calculations on 
Kr atom were performed using the following two basis sets: (i) smaller (20s, 15p, 9d) basis 



set of Koga et al. 



and the (ii) large (32s, 29p, 20<i) universal basis set of Malli et a/.[22| 



Calculated total CPs of Kr, for a.u.< q < 30 a.u., are presented in table HVl which 
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Figure 5: Orbital Compton profiles of Ar, plotted as functions of the momentum transfer q. In 
the decreasing order of the value of J QX \,{q = 0), the profiles correspond to 3s, 3p 3 / 2 /3]?i/2, 2s, 
2p 3 / 2 /2pi/2, and Is orbitals. Note that for all the cases, profiles of p 3 / 2 and pi/ 2 orbitals are 
virtually identical. These profiles were computed using the universal basis set[22j]. 



also contains the nonrelativistic HF profiles computed by Clementi and Roett 



DHF 



profiles calculated by Mendelsohn et a/.[lJ], and the experimental results of Eisenberger and 



Reed 



25]. 



In Figs. [Sandd respectively, our total and orbital CPs of Kr, are plotted as a function of 
the momentum transfer q. Upon comparing our CPs of Kr obtained using two basis sets we 
note that: (i) for small values of q, the values obtained using the smaller basis set of Koga et 
al. 2i| are slightly smaller than the ones obtained using the universal basis set, and (ii) for 
large values of q, the results obtained using the two basis sets are in excellent agreement with 
each other. Next, we compare our calculated CPs with those computed by Mendelsohn et 
a/.fl4| using the numerical orbitals obtained in their DHF calculations. From table IIVI it is 
obvious that, for the all the q values for which Mendelsohn et al. Jj| reported their CPs, our 



22j, are in exact agreement with their results. 

f Mendelsohn et al.[lj|, and our 



profiles obtained using the universal basis set 
As a matter of fact, the agreement between the results o 
results computed using the smaller basis set of Koga et al. 

Upon comparing our results to experimental ones, we see that our universal basis set value 
of J(q = 0) = 7.187, is in excellent agreement with the experimental value of 7.205 25]. 



281 ] . is also excellent. 
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Table IV: Our results on total profiles of Kr computed using the smaller basis set of Koga, Tatewaki 
and Matsuoka (KTM) 
nonrelativistic HF results 



281 ] , an d the universal basis set 



2J] . Re 



241 ] . and the experimental results 



ativistic results of other authors fl^. 



251 ] are also presented for comparison. 



g(a.u.) 


J(q){KTM) a 


J(q)(\Jm) b 


J(g)(DHF) c 


J(q)(RF) d 


J(g)(Exp.) e 


0.0 


7.1788 


7.1871 


7.19 


7.228 


7.205 


0.1 


7.1470. 


7.1548 


7.15 


7.194 


7.152 


0.2 


7.0452 


7.0505 


7.05 


7.085 


7.022 


0.3 


6.8588 


6.8595 




6.888 


6.767 


0.4 


6.5780 


6.5735 


6.57 


6.595 


6.459 


0.5 


6.2087 


6.2010 




6.216 


6.098 


0.6 


5.7744 


5.7670 


5.77 


5.776 


5.701 


0.7 


5.3093 


5.3053 




5.309 


5.289 


0.8 


4.8485 


4.8486 




4.848 


4.880 


0.9 


4.4197 


4.4225 




4.420 


4.491 


1.0 


4.0395 


4.0429 


4.04 


4.039 


4.133 


1.2 


3.4425 


3.4432 




3.441 


3.540 


1.4 


3.0368 


3.0353 




3.037 


3.122 


1.6 


2.7662 


2.7650 




2.769 


2.850 


1.8 


2.5787 


2.5785 




2.583 


2.670 


2.0 


2.4362 


2.4364 


2.44 


2.441 


2.533 


2.5 


2.1425 


2.1428 




2.144 


2.219 


3.0 


1.8571 


1.8572 




1.857 


1.898 


3.5 


1.5784 


1.5782 




1.578 


1.597 


4.0 


1.3257 


1.3255 


1.33 


1.326 


1.338 


5.0 


0.9333 


0.9335 




0.934 


0.937 


6.0 


0.6773 


0.6773 


0.677 


0.678 


0.683 


7.0 


0.5118 


0.5118 




0.512 


0.522 


8.0 


0.4001 


0.4001 




0.400 


0.399 


9.0 


0.3205 


0.3205 




0.319 


0.316 


10.0 


0.2608 


0.2608 


0.261 


0.259 


0.254 


15.0 


0.1062 


0.1062 




0.104 


0.095 


20.0 


0.0506 


0.0506 


re 


0.049 


0.044 


or: n 


n no^7i 


n no^7i 


n no ^7 








Figure 6: DHF Compton profile of Kr, J(q), computed using the universal basis set [22], plotted as 
a function of the momentum transfer q. 



For other values of momentum transfer in the range 0.1 a.u.< q <1.0 a.u., although the 
agreement between our results and the experiments is slightly worse, yet our results are 
closer to the experimental value as compared to the nonrelativistic HF results (3]. For higher 
values of momentum transfer, our DHF results are fairly close to the HF results suggesting 
that in the region of large q, relativistic effects are unimportant. Thus, we conclude that 
from Kr onwards, relativistic effects make their presence felt in the small q region. 

Finally, we investigate the orbital CPs of Kr in Fig. which presents the plots of the 
profiles of outer orbitals starting from 3d 3 / 2 to 4p 3 / 2 . As far as the general trends of the 
orbital profiles are concerned, they are similar to what we observed for the cases of Ne and 
Ar, except for one important aspect. Unlike the Ne and Ar, for Kr for the first time we 
begin to observe the fine structure splitting in the orbital profiles of 4p 3 / 2 and 4pi/ 2 orbitals 
in the low q region, as is obvious from Fig. [7] For example, for q = 0, corresponding 
values are Ji P3/2 = 0.508, and J^ Pl/2 = 0.496, amounting to a difference of ~ 2%. This is in 
complete agreement with our earlier observation that the relativistic effects make significant 
contributions to the CPs of Kr in the small q region. 



17 



0.8 



0.6 




Figure 7: Orbital Compton profiles of Kr for 4s, 4p 3 / 2 , &p\/2, 3s, 3p 3 / 2 /3pi/2) and 3d 5 / 2 /3d3/2 
orbitals in the order of decreasing values at q = 0. For small q values, the differences between the 



4^3/2 and profiles are visible. These profiles were computed using the universal basis set 



22]. 



E. Xe 



In this section, we discuss our results on the relativistic Compton profiles of Xe. The 
DHF calculations on Xe atom were performed using the following two basis sets: (i) smaller 
(22s, 18p, 12d) basis set of Koga et al. and the (ii) large (32s, 29p, 20d) universal basis 



set of Malli et al. 



221 ] . Total CPs of Xe, for selected values of momentum transfer in the 



range a.u.< q < 100 a.u., are presented in table [V] For the sake of comparison, the same 
table also contains DHF, and the nonrelativistic HF, profiles calculated by Mendelsohn et 

n 

al. 14 . Here, we are unable to compare our results with the experiments, because, to the 
best of our knowledge, no experimental measurements of the CPs of Xe exist. 

Additionally, in Figs. [8] and [9] respectively, we present the plots of our total and orbital 
CPs of Xe. Upon comparing our total CPs obtained using the two basis sets we find that, 
as before, they disagree for smaller values of q, with the CPs obtained using the smaller 
basis set [23] being slightly lower than those obtained using the universal basis set 22]. As 
is obvious from table [V], that for q > 1.5 a.u., the two sets of basis functions yield virtually 
identical results. In the same table, when we compare our results to the earlier DHF results of 

n 

Mendelsohn et a/.[14|, we find that for all the q values, the agreement between our universal 
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Table V: Total CPs of Xe computed using the smaller basis set of Koga, Tatewakd and Matsuoka 
(KTM) 2*1. and the universal basis set 



tivistic HF results 



22l |. Relativistic results of other authorsjljl], and nonrela- 



1J| are also presented for comparison. 



g(a.u.) 


J(q)(KTM) a 


J(q)(Vm) b 


J(g)(DHF) c 




0.0 


9.722 


9.737 


9.74 


9.88 


0.1 


9.673 


6.687 


9.69 


9.82 


0.2 


9.515 


9.523 


9.52 


9.65 


0.4 


8.784 


8.775 


8.78 


8.85 


0.6 


7.597 


7.587 


7.59 


7.62 


1.0 


5.448 


5.451 


5.45 


5.46 


1.5 


4.293 


4.292 


4.29 


4.31 


2.0 


3.678 


3.678 


3.68 


3.69 


4.0 


1.707 


1.707 


1.71 


1.72 


6.0 


1.060 


1.061 


1.06 


1.06 


10.0 


0.5150 


0.5150 


0.515 


0.515 


25.0 


0.0660 


0.0662 


0.066 


0.064 


50.0 


0.0088 


0.0088 


0.0088 


0.0076 


100.0 


0.00067 


0.0067 


0.00068 


0.00043 



a our DHF results computed using the basis set of Koga, Tatewaki and Matsuoka 



''our DHF results computed using the universal basis set 



"DHF results of Mendelsohn et al. 
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28]. 



14|] basec 



^Nonrelativistic HF results reported in Ref. 



u 



pon finite- difference calculations 



14] 



basis-set based CPs, and their results, is perfect up to the decimal places reported by them. 
This again points to the correctness of our calculations. 

Upon comparing our DHF results to the nonrelativistic HF results of Mendelsohn et 



al. 14], we find that for smaller values of q, the DHF values of CPs are smaller than the HF 
values, while for large values of q, the trend is just the opposite. 

Finally, upon examining the orbital profiles presented in Fig. [U we observe further 
evidence of the importance of relativistic effects in Xe. As is obvious from the figure, the 
fine-structure splitting between the orbitals profiles of 5p 3 / 2 and 5pi/ 2 orbitals is larger as 
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3, 



Figure 8: DHF Compton profile of Xe, computed using the universal basis set[22j|, and plotted as 
a function of the n 




T 



1 



Figure 9: Orbital Compton profiles of Xe for 5s, 5p 3 / 2 , 5p]/ 2 , and 4d 5/ / 2 /4ci 3 / 2 orbitals in the order 
of decreasing values at q = 0. For small q values, the differences between the 5p 3 in and 5p±/ 2 profiles 



are quite significant. These profiles were computed using the universal basis set 



22] 



compared to 4p 3 / 2 /4pi/ 2 splitting in Kr, and persists for a longer range of q values. For 
smaller values of g, J^p 3/2 {q) > J5 Pl/2 (q), while for large q values, opposite is the case. For 
q = 0, Ji 



5P3/2 



0.592, while J, 



5pi/2 



0.562, which amounts to a difference of 



/o. 
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Table VI: Total CPs of Rn computed using the smaller basis set of Koga, Tatewaki and Matsuoka 



(KTM) 



29], and the universal basis set 221 ] . compared to the earlier calculations of Biggs et al. (if 



<?(a.u.) 


J(q)(KTM) a 


J(q)(XJm) b 


J(<?)(DHF) C 


0.0 


11.8344 


11.8531 


11.9 


0.1 


11.7850 


11.8026 


11.8 


0.2 


11.6176 


11.6306 


11.6 


0.4 


10.8055 


10.7996 


10.8 


0.6 


9.4877 


9.4744 


9.47 


1.0 


7.2130 


7.2130 


7.21 


1.6 


5.8127 


5.8132 


5.81 


2.0 


5.1530 


5.1533 


5.15 


4.0 


2.8379 


2.8381 


2.84 


6.0 


2.0453 


2.0454 


2.05 


10.0 


0.9804 


0.9804 


0.98 


30.0 


0.1083 


0.1083 


0.11 


60.0 


0.0166 


0.0166 


0.017 


100.0 


0.0037 


0.0037 


0.0037 



a our DHF results computed using the basis set of Koga, Tatewaki and Matsuoka 



22] 



29]. 



b our DHF results computed using the universal basis set 

n 

C DHF results of Biggs et a/.[15j] based upon finite- difference calculations 
F. Rn 



As far as atomic Rn is concerned, to the best of our knowledge, no prior experimental 



studies of its Compton profiles exist. However, Biggs et al. 



15|did perform DHF calculations 



of this atom, using a finite difference approach, with which we compare our results later on 



in this section. Our DHF calculations on Rn atom were performec 
basis sets: (i) smaller (25s, 21p, 15c?, 10/) basis set of Koga et al. 



using the following two 
, and the (ii) large (32s, 
29p, 20rf, 15/) universal basis set of Malli et al. 22]. Total CPs of Rn, for selected values of 
momentum transfer in the range a.u.< q < 100 a.u., are presented in table IVTl 

Our results for total and orbital CPs of Rn are plotted in Figs. [U]] and [EH respectively. 
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Figure 10: DHF Compton profile of Rn, computed using the universal basis set|22j, and plotted as 
a function of the momentum transfer q. 

As for other atoms, we find that our total CPs obtained using the two basis sets disa gree 



for smaller values of q, with the CPs obtained using the smaller basis set of Koga et al. 



29] 



22fl. From table 



being slightly smaller than those obtained using the universal basis set 
we deduce that for q > 4.0 a.u., the two sets of basis functions yield virtually identical values 
of CPs. In the same table, when we compare our results to the earlier DHF calculations 
of Biggs et al. 15], we find that for all the q values, the agreement between our universal 



basis-set based CPs, and their results, is perfect up to the decimal places reported by them. 

Of all the rare gas atoms considered so far, on the intuitive grounds we expect the rela- 
tivistic effects to be the strongest in Rn. Indeed, this is what we confirm upon investigating 
the orbital profiles presented in Fig. [HJ As is obvious from the figure, the splitting between 
the orbitals profiles of 6p 3 / 2 and 6px/ 2 orbitals is quite big, and persists for a large range of 
q values. Similar to the case of Xe, here also for smaller values of q, J& P3/2 (q) > Js Pl/2 (q), 
while for large q values, opposite is the case. For q = 0, Je P3/2 = 0.644, while Jq Pi/2 = 0.551, 
amounting to a difference of ~ 15%, which is quite substantial. The fine-structure splitting 
between the profiles of hd$/2 and hd$/2 orbitals although is not quite that large, yet it is 
visible in Fig. CEU At q — 0,Jsd 5/2 = 0.185, and J<od 3/2 = 0.179, leading to a difference of 



~ 3%, which is quite significant for an inner orbital. Thus, we conclude that the relativistic 
effects are quite substantial in case of Rn, and, therefore, it will be useful if experiments are 
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q (a.u.) 



Figure 11: Orbital Compton profiles of Rn for 6s, 6p 3 / 2 , 6pi/2, and 5d5/ 2 /5d 3/ / 2 orbitals in the 
order of decreasing values at q = 0. For small q values, the differences between the 6p 3 / 2 and 6p\/2 
profiles are quite large. Even the splitting of 5d^/2 and 5d 3 / 2 profiles is visible. These profiles were 



computed using the universal basis set 



22] 



performed on this system to ascertain this. 



G. Z dependence of relativistic effects on Compton Profiles 

In earlier sections, while discussing relativistic effects on Compton profiles, we noticed 
that they were most prominent for small momentum transfers. Moreover, one intuitively 
expects the relativistic effects to increase with increasing atomic number Z. In this section 
our aim is to perform a quantitative investigation of relativistic effects on quantum profiles, 
as a function of Z, for both large and small values of momentum transfer. We noticed 
that for small momentum transfers, DHF values of J were smaller than their nonrelativistic 
counterparts, while for large momentum transfer opposite was the case. Therefore, for a 
given value of momentum transfer q, we quantify relativistic effects in terms of | J(DHF) — 
J(HF)|, which is the magnitude of the difference of relativistic DHF value of the Compton 
profile (J(DHF)), and the nonrelativistic HF value of the profile (J(HF)). We obtain J(HF) 
by using a large value of the velocity of light (c = 10 4 a.u.) in the DHF calculations. We 
explore the dependence of this quantity on Z, for two values of momentum transfer, q = 0, 
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Figure 12: Difference between the relativistic (J(DHF)), and the nonrelativistic (J(HF)) Compton 
profiles plotted, on a logarithmic scale, as a function of the atomic number Z. Plots correspond to 
the momentum transfer values q = 0, and q = Z. 

and q = Z a.u., where the latter value clearly belongs to the large momentum transfer 
regime. The values of ln|J(DHF) — J(HF)| as a function of lnZ, are presented in Fig. 
IIH Gl for both these values of momentum transfer. From the figure it obvious that, to a 
very good approximation, the corresponding curves are straight lines, suggesting a power- 
law dependence of the relativistic effects on Z. The slopes of the least-square fit line for 
q = is 2.36 while for q = Z , the slope is 1.35. Of course, these results are based upon 
data points generated by six values of Z (rare gas series), and consequently can only be 
treated as suggestive. But the results suggest: (i) super-linear dependence of the relativistic 
effect on quantum profiles in both momentum transfer regimes, and (ii) stronger influence of 
relativity in the small momentum transfer regime as compared to the large one. Of course, 
this exploration can be refined further by separately investigating the Z dependence of these 
effects on the core and valence profiles. Additionally, this investigation can be extended to 
a larger number of atoms to obtain a larger set of data points. However, these calculations 
are beyond the scope of the present work, and will be presented elsewhere. 
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IV. CONCLUSIONS AND FUTURE DIRECTIONS 



In this paper, we presented an approach aimed at computing the relativistic Compton 
profiles of atoms within the DHF approximation, when the atomic orbitals are represented 
as linear combinations of kinetically-balanced set of Gaussian functions. The approach was 
applied to compute the CPs of rare gas atoms ranging from He to Rn, and results were 
compared to the experimental profiles, and theoretical profiles of other authors, wherever 
such data was available. Additionally, the influence of size and type of basis set was examined 
by performing calculations on each atom with two basis sets: (i ) a well-known smaller basis 
set, and (ii) a large universal basis set proposed by Malli et al. 22]. 

Upon comparing our results with the experiments, we found that for lighter atoms He, 
Ne, and Ar, the agreement was similar to what one obtains from the nonrelativistic HF 
calculations, indicating lack of any significant relativistic effects for these atoms. For Kr, 
we noticed that for smaller momentum transfer values, DHF results were in better agree- 
ment with the experiments, as compared to the HF results. For heavier atoms, Xe and 
Rn, unfortunately no experimental data is available. Yet another quantitative indicator of 
the importance of relativistic effects is the fine-structure splitting of the profiles, i.e., the 
difference in the profiles of ^1/2/^3/2 etc., which will have identical profiles in nonrela- 
tivistic calculations. We found that this splitting becomes larger with the increasing atomic 
number of the atom, thus justifying a relativistic treatment of the problem for heavy atoms. 
Additionally, by comparing our results with the nonrelativistic HF results we found that the 
relativistic effects are most prominent in the region of small momentum transfer, while at 
large momentum transfer, their contribution is much smaller. 

In the literature, we were able to locate prior theoretical calculation of relativistic CPs of 
atoms only from one group, namely the DHF calculations of Mendelsohn et al. lj| and Biggs 
et al. [ll]], performed on Ar, Kr, Xe, and Rn, employing a finite-difference based approach. 
The CPs computed by themQ Q| for these atoms were found to be in perfect agreement 
with our results computed using the universal basis set. This testifies to the correctness of 
our approach, and suggests that by using a large basis set, it is possible to reach the accuracy 



22j, but 



of finite-difference approaches in relativistic calculations, not just on total energies 
also on expectation values. 

Having investigated the influence of the relativistic effects, the next logical step will 
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be to go beyond the mean-field DHF treatment, and incorporate the influence of electron 
correlations on atomic CPs, within a relativistic framework. Such a treatment can be within 
a relatWstic CI f—A o r can a, S o be pe rform e d within a pe rtUrb a t io„- tn eo reti c 
formalism. Work along these lines is currently underway in our group, and the results will 
be submitted for publication in future. 



Appendix A: A DERIVATION OF COMPTON PROFILE MATRIX ELEMENTS 
OVER KINETIC ALLY BALANCED GAUSSIAN BASIS SETS 

During our discussion here, we use the same notations for various quantities as adopted 
in section [TO Our aim here is to evaluate the closed form expressions for the following two 
integrals 

1 r°° 

J t^ = 2J g P9^P)9 L KJ (p)dp (Al) 

J t(v) = 2j q P9 S M9%{p)dp (A2) 

which, as explained in section HD, are needed to compute the orbital (and total) atomic 
CPs when the KBGF based numerical formalism is employed to solve the DHF equations. 
First, we will obtain expressions for g^(p) and g^(p), the radial Fourier transforms of the 
large and small component basis functions g^ir) and gfj(r), respectively, defined as 

9 Ki {p) = ( 27r )3/2 J Q r 9 L nir)j lA { V r)dr (A3) 

roo 

9*i(p) = ,n W2 / r 9i(r)ji B (pr) dr (A4) 



(2tt) 3 / 2 Jo 

where ji A (pr) / ji A (pr) refer to the spherical Bessel functions corresponding to the orbital 
angular momentum Ia/Ib of the large/small component. The spherical Bessel function is 
related to the Bessel function by the well-known relation 

3u(x) = J^Jv+i/2{x), (A5) 
where J u (x) is the Bessel function. 
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1. Derivation for the Large Component 



First , we obtain and expression for g^(p) by performing the integral involved in Eq. 
( 1A3I ). Substituting the expression for g%(r) from Eq. {3} in Eq.jAl), we obtain 



A(p) = 72^71 l°°N^^e-^ 2 j lA (pr)dr 



Kl 



N yn K +l/2) e - ai r J /A+1/2 ( pr ) dr (A6) 



where in the last step, we have used Eq. (IA5I) . Next, on using the relation n K — I a + 1, 
and the definite integral 30| 

x u+1 e~ ar2 J v (J3x) dx = —^-^ e-^/ Aa [Re(a) > 0, Re(u) > 0], (A7) 
the Eq.f lA6l) simplifies to 

9 L M = N -j^m e ~ p2/Aai - ( A8 ) 

On substituting the above result in Eq JAll ). one obtains 

where a i7 = . Next, on making the change of variable t = -p— in the integral above, 

2 

leading to the lower limit q t = we obtain 

(4a jC ^+ 3 / 2 J qt 

leading to the final expression 



T L ' K ( n \ — V kU_\ KjJ V^HjJ / A A -t j-l 

Ji 3 W ~ A {A„,.„,.\U+3/2 I 1 e al > 



4 , {q) Jl^_^_ r(lA + liqtl (A9) 

where T(l^ + 1, qt) is the incomplete gamma function. Since, Ia is a non-negative integer, 
the incomplete gamma function can be easily computed using the series J30]. 

Ia „ m 

T(lA + l,q t ) = (lA)\e^Y,^T- ( A1 0) 
We note that our general result for J^' K (q) in Eq. ( 1A9I ) leads to the same formulas as 



reported by Naon et al.[31j for the atomic CP matrix elements for s- and p-type Gaussian 
orbitals, for the nonrelativistic case. 
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Table VII: Relationship between quantum numbers k, n K , and Ib, for relativistic atomic orbitals. 



K 


n K 


Ib 


~(j + k) 


— K 


j + \ = - K 


U + |) 


K+l 


j - \ = K ~ 1 



2. Derivation for the small component 

Noting that the explicit form of the small component basis function g^(r) (cf. Eq. (j4])) 



is 



<&(r) = iV5JV£ [K + k^-V^ - 2a t r n » +1 e-" 

On substituting the above in Eq. (IA4|) . the Fourier transform of the small component basis 
function becomes 



N S -N L <'>- 

<jU.i>) - ' " 



_ , {n K + K)r^~ l/2) -2a l r^ +z,2) ]e- a ^J lB+l/2 {pr)dr (All) 
VP Jo 1 J 

As before, we seek a relation between n K and Ib, which is summarized in table IVIIL 

Here, the two cases have to be dealt separately since, the relations are different for the two 

possibilities. 



Case (i) K = -(j + 1/2) : 



From table IVIIL it is easy to see that for this case, n K = l B = —k. The integral in 
Eq. llAllI) becomes 



g S M = -^=^ /_ (-2a i )r ( ^ +3/2) e- Q ^ 2 J„ K+1/2 (pr)rfr, 



VP Jo 



-N S -N L - P 

Kl Kl , 



-p 2 /4ai 



(2a i )(^+ 1 /2) 
Case (ii) K= (j + 1/2) : 

For this case, Ib = n n~ 2 = k — 1, which upon substitution in Eq. (|All|) yields 



(A12) 



o N S -N L r°° 

yfP JO 



{2n K - ly^-W - 2a ir ^ +3 ^] e~ a ^ J n ^ 3/2 {pr) dr. (A13) 
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Next we use the result [30I] 



are 



J u ((3x) dx 



(f + 1 + i 



2^+ l «5(^+ 1 )r(z/ + 1) 

for i?e(a) > 0, + 1/) > -1 



* + A*+ 1 « /? 2 
, z/+ 1, - — 

2 ' ' 4a 



(A14) 



where $(0,6,2) is the confluent hypergeometric function, in Eq. (1A13I) . and after some 
simplifications obtain 



tl 



J.) 



(fiK-3/2) 



2 (^-3/2)q;.(^-i/2) ^ k 2 

f l \^f 1 1 P 2 N 

- n K <P n K -\ — , n K , 

V 2/ 1 2 2' 4a,-, 



1 p 2 



2' 4a ?: 



(A15) 

Next, we use the following two identities involving the confluent hypergeometric functions 



and 



a$(a + 1, b, z) = (z + 2a - 6)$(a, 6, z) + (b - a)$(a - 1, b, z), 



$(a, a, = e 2 



to obtain the following simple expression from Eq. (1A15I) 



srJp) = 



S atL V 



-p 2 /4oi 



(A16) 



(A17) 



(A18) 



( 2a .)K+i/2) 

Comparing the results of two cases (1A12I ) and (1A18I ). we find that they only differ by a 
sign, and hence when substituted in the expression for J--' K {q) in Eq.flSS) yield the same 
result 



4^) = 5 / (^)(^)(^)(^) (4 J ai)Hw+1/a ^ /4 ^ ^ 



where a,- 



The above integral can be evaluated in exactly the same way as was 



done before for the large component (cf. IA9I) . to yield the final expression for the Compton 
profile matrix element 



« W 4 (4a,« i )^+ 3 / 2 



r(u + 2, 9t ), 



(A19) 



29 



where q t = y~ , and the incomplete gamma function is defined in Eq. ( 1A10I ). Finally, the 
large and small components of the CP of an orbital can be computed in terms of these 
matrix elements, as 

JnM=Y.C^C^\ (A20) 

JnM=T, C ^ J t- (A21) 



'■J 



It is these formulas derived here which have been numerically implemented in our com- 
puter program COMPTOnQ aimed at calculating relativistic atomic CPs. 
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